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Abstract 

We have carried out numerical simulations of the three-dimensional Ising spin 
glass model with first neighbour Gaussian couplings using three replicas for each 
sample of couplings. We have paid special attention to the measure of two types of 
Binder cumulant that can be constructed from the three possible overlaps between 
the replicas. We obtain new information about the possible phase transition and 
perform an initial analysis of the ultrametricity issue. 



1 Introduction 



The problem of the existence or not of a phase transition in 3-dimensional spin glasses 
is still a subject of controversy 0. In the last year a lot of work has appeared with the 
purpose of clarifying the existence or not of the phase transition, and in the first case, 
characterize the critical exponents 0, |3|, |Q. Moreover recent numerical work has shown 
that the low temperature phase is as predicted by Mean Field ||. Also recent analytical 
work by F. Guerra || has clarified the meaning of some formulae found in the framework 
of the Mean Field Theory. 

The most relevant observable in the study of these systems is the overlap between 
replicas, i.e. the order parameter of the system. Up to now, all simulations have been 
carried out with only two replicas for each sample of couplings. The principal novelty we 
introduce in this paper is the use of three replicas. 

A way to study a possible phase transition is to use the Binder cumulant. This quantity 
clearly marks the change to a Gaussian situation from the non Gaussian one. At present 
this approach has not been very successful for spin glasses because at low temperatures 
the results turn out to be inconclusive. One needs a lot of statistics and large lattices to 
obtain the small deviation between curves corresponding to different sizes. 

In this paper we study two types of Binder cumulant constructed from the three different 
overlaps we can measure between the three replicas. The objective is to obtain a clearer 
signal than in the case of two replica cumulant. From one cumulant, which shows a signal 
very similar to the two replica ones, we have estimated the critical temperature and the 
exponent v on the paramagnetic side, obtaining a good agreement with the previous quoted 
values @ . The other cumulant presents new interesting features and supplies us a method 
of studying the possible ultrametricity of the spin glass phase as well as a determination 
of the critical temperature. 

2 Model, simulation and observables 

The model we have studied is the 3 — d Ising spin glass with nearest neighbor couplings 
distributed Gaussianly around zero. The Hamiltonian is 

H = - JijOiOj ■ (1) 

<i,j> 

As usual < i,j > denotes nearest neighbor pairs. The lattice sizes we simulated are 
L = 4,6,8,10. We have used the simulated tempering method [0, §]. The range of 
temperatures studied has been [0.7, 1.3] (for L = 4 we have also made a run in the range 
[0.3,1.3]) in steps of 0.05. We have simulated 4096 different samples of couplings for 
L = 4, 6 and 2048 samples for L = 8, 10, with three replicas for each sample^. 

The number of sweeps slightly depends on the lattice size, but it is typically one million 
for measuring, after the order of half a million iterations to estimate the free energy. The 

1 Let be a, r and \i the three replicas that we will simulate in parallel with the same disorder. Hence, 
we can define three different overlaps that we will denote {912,923,913} or {q,q', q"} indiscriminately in 
the rest of the paper. 
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calculations have been carried out on a tower of APE100 |J with a real performance, for 
this problem, of 5 Gigaflops for a total time of three weeks. The errors, sample to sample, 
have been computed with the jack-knife method. 

As a check of the thermalization procedure we have monitored the symmetry of the 
distribution of the overlaps (i.e. P(q) = P(—qj). We have also checked the two following 
relations |§: 

Wf = \W) 2 + \W> (2) 
1- 



<?y 2 > = 2<4 2 r+2<9 4 > (3) 

As usual, we denote thermal averages by ((••)) an d disorder averages by (••) The first 
formula (fj) was pointed out in a previous numerical analysis ||, and have been rigorously 
demonstrated, along with (|j), by F. Guerra ||. Both relations are well satisfied by our 
data for every value of L and T. Their validity has been proved in the infinite volume 
limit, however it is possible that the deviations are small also for not too large volumes, as 
those of the present paper. 

We have focused our attention on the following Binder cumulants constructed from the 
three overlaps measured for each sample of couplings: 



R = (1 212913423 1) R/ = (212913923) m 
^qqq — 3 / 2 ° 111 ~ 3/2 V / 

(q 2 ) ' (qT 



and 



R = (QgigJ - kisl) 2 ) R , = ((912-913 sign(g 23 )) 2 ) 

D q-q — -, ff q-q — 7 \°) 

(923 2 ) (923 2 ) 

where ^23 is the largest one (in absolute value). These definitions must follow the finite 
size scaling relation 

B # = f # (L 1/ »(T-T c )) (6) 

where we have used the symbol # to denote either cumulant. In the following sections we 
will study, in detail, these observables. 



^ Bqqg cUld B gqg 

In Fig. 1 the values obtained from the simulations for B qqq are shown (lower curves). 
We observe that at high temperature the values for the different lattice sizes are widely 
separated (the larger L the smaller B qqq ). However below T ^ 1.1 the curves intermingle. 
In order to distinguish them clearly it would be necessary to reduce considerably the errors. 
The observable B' qqq presents a similar behavior. In fact, the product 912913923 is generally 
positive (because it is favored combinatorially even if one has three completely independent 
configurations) . 

In the same figure, we show the values we would expect if the three overlaps satisfied 
ultrametricity (upper curves). We have calculated these curves in the following way. The 
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Figure 1: B qqq against T: Numerical data (lower curves) and supposing ultrametricity (as 
explained in the text) (upper curves). 
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probability distribution of the overlaps for each temperature and lattice size, P(q), is 
obtained directly from the simulations. From it we extract the function x(q) defined as 

x{q) = f 9 dq'P( q >), (7) 



o 

and then we obtain the inverse q(x), which is shown in Fig. 2 for different temperatures 
at L = 10. Supposing ultrametricityQ |1(J 

Ps(q, q', q") = \p{q)x{q)8{q - q')8{q - q") + 

- {P{q)P{q')0{q' — q")5{q" — q) + two permutations) (8) 

the cumulant will take the value 

1 /3 r 1 r 1 1 r 1 \ 

B ig<i = v~, ^3/2 o / dx( i 2 ( x ) / d y<i(y) + o / dxx( i 3 ( x ) ■ ( 9 ) 

We can observe that the curves extracted directly from the simulations approach those 
obtained supposing ultrametricity when the temperature decreases. This is in agreement 
with the fact that ultrametricity could be expected at low temperatures. However it is not 
demonstration because practically any reasonable type of relation (including independence) 
between the q's would produce curves approaching the value 1 when the temperature goes 
to zero. On the other hand, the curves from the simulations are under the value 1 (the 
points for L = 4, T = 0.3 — 0.4 are over 1 but contain this value inside their error 
bars) while the lines bases on ultrametricity always remain greater than 1. These results 
do not rule out the possibility that the overlaps are as independent as possible for any 
temperature (even when the three replicas are totally independent, the probability of the 
three overlaps has a constraint given by a simple combinatorial problem) ; when decreasing 
the temperature the behavior of the P(q) favors the B qqq approaching 1. 

The equation @ can also be used to see what would be obtained if we assume for q(x) 
a shape that could be expected in the thermodynamic limit, that is to say 

q ( x ) = Slo(i}\ 0<x<m, (1Q) 
[ q , otherwise. 

We would obtain 
B, 



qqq ~ H 2r 
V 2r+l 



rn 



)3/2 




2r + 1 



In Fig. 2 we observe that the function q(x) obtained from the simulations is rather 
similar to the supposed here with m = 1 and r descending from approximately 1 for 
high temperature towards as the temperature diminishes. As a test, we can compare 
qualitatively the curves shown in Fig. 1 in the case of ultrametricity with a function of 
this type, and it is seen that in fact the behavior is similar. 

2 In the appendix we will show that if ultrametricity, convexity and positivity hold then the functional 
form of Ps{q, q' , q") must be that of Mean Field for a generic P(q). 
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Figure 2: q(x) for 13 different temperatures. 
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Figure 3: Probability to have a b' value for the four different size and T = 0.8. 



We have estimated the value of the exponent v by fitting the derivative of B qqq in the 
high temperature region as 



f(g,L)- dB ™ 



dT 



= aL 1/u . (12) 

To: B qqq (T )=g 



We have used polynomials of different orders to fit B qqq (T) and calculated the derivatives 
at several values of B qqq . Let us remark that these derivatives are taken at different values 
of T but at fixed B qqq . All the results obtained for u, for different values of g near the 
critical one g c = B qqq (T c ), are compatible within the errors, which turn out to be large. 
Finally we give an estimation of v = 1.5(3), which is in good agreement with the value 
v = 1.7(3) reported in reference for the ±J spin glass. 

As we have remarked above, the product gi2Q , i3?23 is mainly positive, but not always. 
In Fig. 3 we show the probability distribution of 

qqq ~ ate 2 2 +^+^ 3 )) 3/2 { ] 
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Figure 4: Scaling, as a function of size, of the area of the negative part (b' qqq < 0) in figure 
3. 



for the different L values at T = 0.8. Fig. 4 attempt to show how the area of the negative 
part behaves with L, but it is not clear if it will go to zero with increasing L or not. If 
mean field holds the negative area must go to zero. 



4 B q _ q and B' q _ q 

As can be seen in Fig. 5, B q _ q has a much clearer signal than B qqq . We observe a crossing 
of the curves for different L's in the region T ~ 1. At low temperatures the values for 
L = 10 mask an inversion of the curves (i.e. the larger L the smaller B q _ q while the high 
temperature order is the other way round) which is very clear for L = 4,6,8, as shown 
in the expanded Fig. 5(6). From the crossing of the curves for the larger lattices we can 
estimate a critical temperature of approximately T c = 1.02(5). B' q - q also presents similar 
behavior. 

These observables, B q _ q and B' q _ q , would be zero if ultrametricity were exactly verified. 
In any case, we expect violations of ultrametricity due to the finite size of our lattice as 
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Figure 5: B q _ q versus T for the four different sizes (upper). In the lower part we plot only 
L = 4,6,8 and T G [0.7,1.3]. 
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happens in the SK model [OH- The same considerations as in the preceding section can be 
done regarding to the approach of the B q _ q curves to zero when the temperature decreases. 
A trial of estimating v like above gives results that move systematically with the value 



B q _ q fixed and then it has been discarded. 



We have made a study of the probability distribution of 

^ (kisi-kalf (14) 

Q23 2 

(where #23 is the biggest one) calculated at every single iteration, at temperature T = 0.8. 

In Fig. 6 we show the logarithm of the probability of having a value of b larger than a 
certain 6 versus 6q. Fig. 7 has the same variable on the vertical axis but the logarithm of 
60 on the horizontal one. The curves for the different values of L are hardly distinguishable. 
Expanding the image in the region of small 60 (figures 6(6), 7(6)), we see that the upper 
curves are those of smaller L corresponding to bigger B q _ q , but the lowering with L is 
slow and it is difficult to guess any asymptotic behavior. Selecting for instance a value 
of 60 = 0.05 and observing how this probability decreases with L we find approximately 
a power law with a small exponent of 0.13(2). In the region of large 6 , the upper curves 
correspond to larger L (in order to maintain the total probability normalized). For all the 
values of L, we observe a first region, for small 60, where the decrease of probability can be 
approximated by a power behavior of the form 6o~ a with a = 0.71(2). There is a central 
plateau where it decreases exponentially as e~ ab ° with a ranging from 5.8(1) for L = 4 to 
5.3(1) for L = 10, while it goes to zero faster when 6 approaches the geometrical limit of 
1. 

The observable 

y _ (gi2 - gig sign(g 23 )) 2 

<?23 2 

shows a similar behavior. 



5 Conclusions 

We have reported in this paper a study of the three dimensional Ising spin glass model 
using three replicas. In particular we have studied different version of the Binder cumulant 
with three replicas. We have estimated both the transition points as well as the critical 
exponent v being in a good agreement with previous reported values 0. 

Using one of the Binder cumulants we have begun a preliminary study of ultrametric- 
ity We have found that the 6 probability distribution goes to zero follows a power law. 
This results are not conclusive (the exponent is small) and we will need simulate large 
lattices using new thermalization methods (for instance parallel tempering 0]) in order to 
distinguish between lack of ultrametricity and possible violations due to finite size effects. 
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Figure 6: Plot of the accumulated probability for all the points (upper) and only for a 
window near the origin (lower). 
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Figure 7: The same of Figure 6 but in a double log scale. 
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7 Appendix 

In this appendix we will show that if ultrametricity holds as well as convexity and positivity 
then the three replica probability must be the Mean Field one. 

The most general formula, using the invariance of the Hamiltonian under the exchange 
of replicas, for the probability of three replicas assuming ultrametricity is: 

P 3 (q, q', q") = A(q)S(q - q')8{q - q") + B(q, q')6{q - q')5{q' - q") + 

B(q', q")6{q> - q")8{q" - q) + B(q", q)B{q" - q)8(q - g'), (16) 

with A(q) and B(q,q') satisfying normalization condition and B(x,y) = B(y,x). Integrat- 
ing q" we obtain the probability to have two replicas 



P2(q,q') 



POO 

A(q)+ / dq"B(q',q" 
J a 



5(q-q') + B(q,q'), (17) 



and finally the probability of one replica {P(q) that we will denote , in this appendix, as 
Fi(<z))is 

roo roo 

P l (q) = A(q)+ dq'B(q,q')+ dq'B(q,q'). (18) 

Jq J — oo 

Now we impose the relation (demonstrated by F. Guerra using only positivity and 
convexity) 

/OO 1 o 1 

dqdq>P 2 (q,q')q 2 q' 2 = -(q*} 2 + -(q*) (19) 
-oo Z Z 

where _ 

roo 

(q-) = / dqP 1 (q)q n . 
We finally obtain the following equations: 

A(q)= [ q dq'B(q,q') (20) 

J — oo 

B(q, q') = 2 QT dq"B(q, g")) (jT dq"B(q', g")) (21) 
Joining equations ([18|) and (p0|) 

/•oo 

P 1 (q) = 2 dq'B(q,q') (22) 



and then we recover the Mean Field formulas, i.e. 

A(q) = l -x{q)PMi (23) 

B(q,q') = ±P 1 (q)P 1 (q') (24) 

with a free P\(q). Obviously this development does not imply any particular functional 
form for Pi{q). 
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